Quadrupole oscillation in a dipolar Fermi gas: hydrodynamic vs collisionless regime 
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The surface quadrupole mode of an harmonically trapped dipolar Fermi gas is studied in both the 
hydrodynamic and collisionless regimes. The anisotropy and long range effects of the dipolar force 
as well as the role of the trapping geometry are explicitly investigated. In the hydrodynamic regime 
the frequency is always slightly smaller than the \/2uj± value holding for gases interacting with 
contact interactions. In the collisionless regime the frequency can be either pretty smaller or larger 
than the non- interacting value 2ijj±, depending on the cloud aspect ratio. Our results suggest that 
the frequency of the surface quadrupole oscillation can provide a useful test for studying, at very 
low temperatures, the transition between the normal and the superfluid phase and, in the normal 
phase at higher temperatures, the crossover between the coUisional and collisionless regimes. The 
consequences of the anisotropy of the dipolar force on the virial theorem are also discussed. 
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I. INTRODUCTION 

Soon after the experimental realization of Bose- 
Einstein condensation and Fermi superfluidity in atomic 
gases, much attention has been devoted to the study of 
collective oscillations (for reviews see, for example, [H,!!]). 
Collective oscillations actually provide unique and valu- 
able information on the macroscopic and dynamic behav- 
ior of these quantum gases, being sensitive to both the 
equation of state, the effects of superfluidity and the role 
of collisions. While compressional modes are directly sen- 
sitive to the equation of state and have been employed to 
test its fine details with high precision (see, for example, 
the case of the Fermi gas along the BEC-BCS crossover 
0), surface oscillations are more suited to understand 
whether the gas is in the hydrodynamic or collisionless 
regime. This is especially important in view of the possi- 
bility of identifying the superfluid phase of the gas, where 
the dynamics of the low frequency modes is always gov- 
erned by the equations of irrotational hydrodynamics. 

The recent experimental availability of atomic gases in- 
teracting with magnetic dipolar forces Q and the very re- 
cent progress in the realization of ultracold hcteronuclear 
molecules interacting with electric dipolar forces Q have 
already stimulated a significant number of theoretical pa- 
pers (for a review see, for example, @). Special attention 
has been also focused on the study of the collective oscil- 
lations of both Bose la and Fermi H-ini dipolar gases, 
inclu ding the study of long range effects in bi-layer sys- 
tems [l^ • A major focus of most of these papers was the 
problem of the instability caused by the attractive part 
of the dipolar force which shows up in the appearence 
of an imaginary value of the collective frequency of com- 
pressional modes. In this paper we will instead focus on 
the main features of the surface oscillations with the aim 
to explore the consequences of the anisotropic and long 
range nature of the dipolar force in both the hydrody- 



namic and collisionless regimes. In this work the dipolar 
effects are taken into account within mean-field theory. 

For harmonically trapped gases interacting with zero 
range forces, which is the case of atomic gases with- 
out dipolar interactions, the frequency of the surface 
quadrupole mode, in the hydrodynamic regime, is given 
by the universal value \^uj± [13| where uj± is the radial 
oscillator frequency in the x — y plane where we excite 
the collective oscillation. This value differs from the non- 
interacting value 2uj± of the ideal gas since in the hydro- 
dynamic regime the kinetic energy does not contribute to 
the restoring force of the quadrupole oscillation which is 
then entirely provided by the external harmonic poten- 
tial. In the collisionless regime the kinetic energy does 
instead contribute to the restoring force and the value 
of the frequency of the quadrupole oscillation depends 
on the explicit balance between the kinetic and exter- 
nal potential contributions. This balance is affected by 
the presence of two-body interactions and is usually ac- 
counted for by the virial relationship [l^ . The anistropic 
and long range nature of the dipolar force affects the 
above scenario both in the hydrodynamic and in the col- 
lisionless regimes. We will mainly focus on the case of 
Fermi dipolar gases where the quest of superfluidity is 
more challenging for single spin species gases and the 
competition between the hydrodynamic and collisionless 
regimes is relevant also at very low temperatures. The 
paper is organized as follows. In Sect. II we discuss the 
equilibrium properties of the gas and the anisotropy ef- 
fects of the dipolar interaction on the density profile. We 
furthermore develop a variational approach, based on a 
scaling ansatz, to derive the equations of motion charac- 
terizing the quadrupole oscillation in both the hydrody- 
namic and collisionless regimes. In Sect. Ill we discuss 
the resulting effects on the collective frequencies and the 
temperature conditions for being in the collisionless or 
hydrodynamic regime. In Sect. IV we provide the ex- 
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plicit derivation of tfie virial tfieorem and we point out 
its relevance for experiments where one suddenly switches 
off the dipolar interaction. In Sect. V we draw our con- 
clusions, while in the Appendix A we present the results 
for a few relevant dimensionlcss functions whose know- 
eldge is needed for the actual calculation of the collective 
frequencies. 

II. TRAPPED DIPOLAR FERMI GASES 

We consider a Fermi gas of dipolar atoms or molecules, 
confined in an axially symmetric harmonic trap 

VU^) = ]^m{ulrl+ulz^), (1) 

where m is the mass of a single fermion, uj±_ and lo^ are 
the radial and axial trapping frequencies, respectively, 
and r^^ = -\~ . If the dipoles are polarized along z, 
the microscopic dipolar interaction can be written as 

^^'i(r-r') = ^^^(l-3cos2 0) , (2) 

with dp = p^/47reo for electric dipoles (p is the electric 
dipole moment) = hq^'^/Att for magnetic dipoles (/i 
is the magnetic dipole moment) The angle 6 is the angle 
that the polarization axis forms with the relative distance 
between two dipoles. 



di it was shown that both the anisotropy of the 
kinetic energy and the exchange effect in the dipole en- 
ergy can be relevant for establishing the exact conditions 
of stability of the gas, but are negligible if one considers 
configurations sufficiently far from the instability point. 
Notice the absence, in Eq.(j31), of the contact term fixed 
by the s-wave scattering length which typically enters the 
energy functional of interacting Bose gases. This term is 
absent in the single species Fermi gas considered in the 
present paper. 

The equilibrium density profile of the trapped gas is 
obtained by imposing minimization of the energy E with 
the proper normalization constraint, which is obtained 
adding the grand canonical term — /i J drn(r) to the en- 
ergy functional, with fj, the chemical potential. In general 
the equilibrium profile will be a function of the quantity 
(p— Vho(r))- In the following, for simplicity, we will make 
the further assumption that the density profile takes the 
form 

where N is the number of fermions and the radii R± 
and Rz are variational parameters. The density Eq.([7]) 
would correspond to the exact Thomas-Fermi solution 
of a Fermi gas in the absence of the dipolar interaction 
term. Using this ansatz, the dipole energy ii^Hip can be 
evaulated analytically and takes the form [3, : 



A. Equilibrium properties 

The ground state properties of the system can be cal- 
culated employing a variational approach starting from 
the following energy functional, based on a mean field 
picture: 



where 



E = = Skin + Sho + Edip 



5 2m 



(3) 



(4) 

-Bho = j drn{r)Vi,^ (5) 

-Bdip = j dYidr2n{Yi)n{r2)vdd{ri - Y2) (6) 

are, respectively, the kinetic, oscillator and dipolar en- 
ergies. In writing the above equations we have used 
a Thomas-Fermi description for the kinetic energy, un- 
der the assumption of isotropy of the momentum dis- 
tribution. This assumption, which applies to three- 
dimensional configurations, neglects possible anisotropy 
effects caused by the dipole interaction in momentum 
space. In the expression for the dipole energy we have 
further ignored the occurrence of exchange terms. In 



£^dip = --^Nepeddlii^)- 



(8) 



In the previous expression ep = h'^kp/2m is the Fermi 
energy, kp being the Fermi momentum defined by the 
normalization condition R^Rzkp = 48N and we have 
introduced the dimensionlcss parameter 



£dd 



d^ikpf 
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(9) 



The dimensionlcss function /(k) is defined in the Ap- 
pendix, with k = R±/Rz the aspect ratio of the trapped 
cloud Notice that in general k differs from the as- 
pect ratio A = ujz/u}± of the trapping potential as a con- 
sequence of the anisotropy of the dipolar interaction. It 
is also interesting to notice that the term in front of /(k) 
in the dipolar energy Eq.® scales like 1/RxRyRz, which 
is the same scaling dependence of a zero-range contact 
force. Deviations from the behaviour of a contact force, 
caused by the anisotropy and finite-range nature of the 
dipolar force, are then entirely accounted for by the func- 
tion /(k). 

The equilibrium configuration is obtained by imposing 
the stationarity of the energy functional Eq.Q with re- 
spect to variations of the Thomas-Fermi radii. One then 
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finds the following equations: 




(10) 
(11) 



where k = y^Ti/moji, i =-L,z is the liarnionic oscillator 
length in the i-th direction. Combining equations (|10p 
and (jlip one finds a transcendental equation for k as a 
function of A and Edd, 



K 

A2 



1 - K 



1 



2edd - 1 



edd ~ 1 (12) 



Note that this expression is the same as the one found 
for bosons [20| . provided the parameter Edd is properly 
defined. 

It is useful to express the previous quantities in terms 
of experimentally controllable parameters. We define the 
non-interacting Fermi energy = h{6Nu}'\^ujzY/^ , the 
non-interacting Fermi momentum kp = y^2me^/'h and 
the corresponding dimensionless parameter 



d^(fc^) 



\3 
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(13) 



In Figure [T] we show the relationship between the rele- 
vant interaction parameter Edd and e'^^. Since £dd/£dd ~ 
kp/kp, the difference between the two quantities is due 
to the departure of kp from the non-interacting value kp, 
caused by the dipolar interaction. 

In Figure[2]we report instead the dependence of k/A on 
Sdd for different values of A. This figure explicitly reveals 
how the dipolar interaction modifies the aspect ratio of 
the trapped gas with respect to the trap value A. For 
zero dipolar interactions we recover the non-interacting 
limit, for which k = A. Note also that, as expected, 
the magnetostriction (electrostriction) effect is larger for 



A 1, while it is much smaller for very small or large 
trapping aspect ratios. The dashed lines in the figure 
mark the regions for which the solution, inclu ding the 
exch ang e interactions, would predict instability 15[ (see 
also [ll[ [ll] ) and we will keep this notation all over the 
work. 



B. Variational approach to the m = 2 surface 
quadrupole oscillation 

In the following we study the effect of the long-range 
dipolar interaction on the frequency of the transverse 
(to = 2) quadrupole mode assuming the gas is either 
in the collisionless (CL) or in the hydrodynamic (HD) 
regime. The results for the latter are always attained in 
the superfluid phase, whose description is based on the 
equations of irrotational hydrodynamics. 

In order calculate the collective mode frequencies we 
implement a variational approach based on the action 



S 



(14) 



from which we can derive the equations of motion for con- 
figurations perturbed from equilibrium and hence evolv- 
ing in time. The m ~ 2 quadrupole oscillation can be 
described employing a variational approach based on the 
scaling transformation 



)T,kii^kPk-ykPl)+hc) i^j,^ 





(15) 

where the two time dependent parameters a{t) and b(t) 
account, respectively, for the presence of current terms 
and changes in the density distribution with respect to 
the unperturbed configuration |^')q. Here is the i — th 
component of the momentum and k = 1,...,N is the 
particle index. The density distribution and the velocity 
field associated with the wave function I^P) are given by 



nix,y,z) = no{e' 



'''(*)a;,e''«y,z) 



(16) 




A,<1 



X>1 



FIG. 1: Density-dependent dipolar parameter edd as a func- 
tion of the externally controllable dipolar paramenter eJJ^j. 




FIG. 2: Gas aspect ratio normalized by the trap aspect ra- 
tio, n/X, as a function of Edd- Left panel: for a cigar-shaped 
trapping potential; Right panel: for a pancake-shaped trap- 
ping potential. Dashed line marks instability as predicted by 
including exchange interactions (see for instance [IF 
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and 



m 



(17) 



where no is the equilibrium density profile. 

The first order variation with respect to a and 6 has 
to be zero, while from the second order variation we find 
the equations of motion for the parameters a and h in 
the linear regime from which we calculate the collective 
frequencies. The term associated with the time derivative 
in (fT4|) gives: 



-ih 



d_ 

"dt 



* ) = 2hahN{r 



(18) 



where the radial square radius (r^) should be calculated 
at equilibrium. 

The hamiltonian term in the action integral (|14p is just 
the energy functional Eq. ([S]). However, out of equilib- 
rium the expression for the kinetic energy differs from the 
Thomas- Fermi expression Eq. ([4]). First of all the pres- 
ence of the velocity field Eq. (fT7| , caused by the scaling 
transformation, gives rise to a "classical" contribution to 
the kinetic energy of the form 



(5-Bkin(ck) 



(irv^n(r) 



2r 



Na^{r\). (19) 



Second one should distinguish whether the system is 
in the hydrodynamic regime, where local equilibrium 
is restored during the oscillation, or in the coUisionless 
regime. In the first case one can still use the Thomas- 
Fermi expression Eq. Since this expression includes 
a local dependence on the density, the density change 
Eq. (jl6p caused by the quadrupole scaling transforma- 
tion does not modify the value of the kinetic energy which 
is consequently affected only by the velocity field contri- 
bution Eq. ([T9|) . In the coUisionless case the kinetic 
energy gets an additional contribution due to the static 
part of the transformation (term in h). The latter causes 
a deformation of the Fermi sphere in momentum space 
giving rise to an elastic type contribution, a typical fea- 
ture of Fermi systems [l3| • The total contribution to the 
changes in the kinetic energy term takes in this case the 
form 



5E, 



(CL) 
kin 



5Ei,Ua)+2b'Ei, 



(20) 



where i?kin,_L is the radial contribution to the ground 
state kinetic energy. 

The other terms in the energy functional (oscillator 
and dipolar energy) depend only on the local density dis- 
tribution and are consequently affected by the scaling 
transformation through the term in h{t): 



8E^,{b)^mulNb'{r\), 
5Edip{h) = --b^Nepeddgin). 



(21) 
(22) 



In order to derive the dipolar contribution Eq. ([22)) we 
have expanded the dipolar energy Eq. ([5]) to second or- 
der in b, using the scaling transformation — > e^Rx and 
Ry — >■ e~^Ry and have introduced the function ^(k) pro- 
portional to the second derivative of £'dip(6) with respect 
to b (see the Appendix A). 

Using Eqs. (|TH1) and (IT11)-(111]) the action in the hydro- 
dynamic regime reads 



5( 



HD) 



+Jdt 

+ 2habN(r^: 



2h' 



Na^{r\) + muj\Nb'^{r\) + 



b^NepEMgin) 



(23) 



where S'^^^ is the action calculated for the ground state. 
In the coUisionless regime one should add the term 
/ dt2b'^EMn,^ to the right hand side of Eq. 

Imposing that the action be stationary with respect to 
variations of a{t) and b{t) the equations of motion can 
be written as 



2h 

2-±-Na{rl) - 2hbN{rl) = 



(24) 



2mu}lNb{rl) + 2haN{r\) - -bNepeddgiK) = (25) 

in the hydrodynamic regime while, in the coUisionless 
regime, they read 



9h 

2 Na{rl) ~ 2hbN(rl) = 

TO 



(26) 



46^;kin,± + 2mujlNb{r'i) + 2naN{r']_) - -bNepEddgin) = 

(27) 

The above equations admit harmonic solutions. The fre- 
quency of the oscillation in the hydrodynamic regime is 
given by 



,(HD) 



= \/2( 



1-1 
4 



1 + £dd 



3 ff/W 
2 



1 



g{i^) 



1/2 



(28) 

which reduces to the usual hydrodynamic result \f2uj\^ 
by setting Edd equal to zero (notice, however, that in a 
single species Fermi gas the Edd = limit is incompatible 
with the hydrodynamic regime which always requires the 
presence of interaction terms). The same expression Eq. 
(j28p for the quadrupole frequency in the hydrodynamic 
regime holds also for bosons, where Edd = m(P /Sh^ag 
with Qs the s-wave scattering length. 

In order to provide a compact expression for the 
quadrupole frequency in the coUisionless regime it is use- 
ful to employ the radial virial relationship (see derivation 
in Sect. pV)) ) 



2E, 



'kin,± 



2E\ 



ho,_L 



2E, 



1 



dip 



NepeddKn) = (29) 
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to express the kinetic energy contribution entering 
the equation of motion (P?]) in terms of the oscilla- 
tor and dipolar energy contributions. Here i^ho.i = 
muj']_Nb{r\)/2 and -Bdip is given by Eq. (P. The func- 
tion h[hi) = —Kdf{K)/dK is proportional to the first 
derivative of the function /(k) entering Eq. ([8|) and is 
explicitly calculated in the Appendix A. Using the virial 
relationship (P^)) one can derive the following useful ex- 
pression for the quadrupolc frequency in the coUisionless 
regime: 



X=0, 



^=4 



(CL) 



1+1. 



-dd 



1 + £dd 



\2 



X ( /(k) + -h{K) - -g{n) 



- 1 

1/2 



(30) 



showing that, as Sdd ^ 0, the non-interacting limit 2uj±^ 
is recovered. 

Eqs. (pS)) and (PU]) represent the main result of the 
present paper. They give explicit predictions for the 
quadrupolc frequencies in the hydrodynamic and coUi- 
sionless regimes in terms of the relevant dimensionless 
parameters Edd and k. These parameters are directly re- 
lated to the experimentally controllable parameters ej]^ 
and A, as shown in Figs. [T]and[2] 



III. QUADRUPOLE MODE IN THE 
HYDRODYNAMIC AND COLLISIONLESS 
REGIMES 

We are interested in comparing the frequency of the ra- 
dial quadrupolc mode in the hydrodynamic and collision- 
less regimes. Let us first analyze Eqs. (pS)) and ([50)) sep- 
arately. In the hydrodynamic regime the quadrupolc fre- 
quency is shifted with respect to the usual hydrodynamic 
result ^/2ijj±, holding for gases interacting with contact 
forces. The deviations are caused by the anisotropic and 
finite range character of the dipolar interaction. This 
shift vanishes for k — > and k — > oo, since in these two 
limits the system is so deformed that the dipole force ef- 
fectively behaves like a contact force. This is reflected in 
the flat behavior of /(k) for very pancake and very cigar- 
shaped traps, as well as in the fact that g{Hi) tends to zero 
in these limits (see the Appendix A). The top panels of 
Fig. [3] show the behavior of the quadrupolc frequency in 
the HD regime as a function of the dipolar parameter Sdd 
for different trap aspect ratios A. The dashed lines mark 
the values of Sdd for which the Fermi gas is no longer sta- 
ble, according to [§] . The figure shows that the frequency 
is always smaller than the ^/2u)±_ value and that the shift 
is largest for intermediate trap asp ect ratios. However, 
as already pointed out in Ref. the shift of the sur- 
face quadrupolc oscillation is in general small compared 
to the case of comprcssional modes which are much more 
sensitive to the equation of state. 




FIG. 3: Top panels: Quadrupole frequency as a function 
of the dipolar strength Edd, both for the CL and the HD 
regimes and for different trap aspect ratios. Bottom panel: 
Quadrupole frequency as a function of A for e^^^^^^^ = 
0.5. Dashed line marks instability as predicted by including 
exchange interactions (see for instance [3). 



In the coUisionless regime the situation is very differ- 
ent. In fact the frequency shift with respect to the ideal 
gas value 2u)^ depends strongly on the geometry of the 
trap and on the strength of the interaction, as can be seen 
from the top panels of Fig. [3l In general we can say that 
for cigar-shaped clouds, the shift is positive, while for 
pancake-shaped clouds the shift is always negative. For 
A ^ 2 wc find that the shift is very small for all values 
of Edd- In contrast to the hydrodynamic regime, in the 
coUisionless regime the frequency shift of the quadrupole 
mode docs not vanish in the limit of large shape deforma- 
tions K (pancake geometry), but instead tends to a value 
that depends only on Edd'- 



,(CL) 



(k — >■ oo) = 2a; J 



1 - 



£dd 



1 + 2Edd 



1/2 



(31) 



For very large Edd Eq. pip approaches the hydrodynamic 
value v^wj^. This is not a surprise since in this limit the 
kinetic energy contribution to the energy of the dipolar 
gas is strongly suppressed, the equilibrium profile being 
the result of the competition between the harmonic po- 
tential and the dipolar force. The large Edd limit is remi- 
niscent of the zero sound scenario in Fermi liquids where, 
for large values of the interaction coupling constant, the 
zero sound velocity approaches the value of first sound. 

The bottom panel of Fig. |3] shows the dependence of 
the quadrupole frequency on the trap aspect ratio, for a 
given value of the dimensionless parameter 



9l2el/6 /^3, , \ 1/2 



34357r2 



(32) 



We have identified this parameter since, for a given value 
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of Lij±, it is independent of the geometry of the gas (A 
or k) or its density. The value of the numerical factor 
in front of (P is 0.1973. To give an order of magnitude, 
using the experimental parameters of Ref. [26j . that is 
N = 2200 ''OK^^Rb molecules in a trap with wj_ = 36 x 2it 
s"^, the condition e^j^A^^/^ = 0.5 corresponds to an elec- 
tric dipole moment of 0.74 Debyc. We see from the figure 
that the frequencies in the hydrodynamic and collision- 
less regimes are sensibly different. They can be used 
therefore to study the transition (or crossover) between 
the two regimes. It can take place in two different scenar- 
ios: when the gas enters the superfluid phase, or when 
the gas goes from the coUisionless regime to the collisional 
regime, where hydrodynamic equations are imposed by a 
high collision rate. 



For very pancake systems, superfluidity cannot take 
place because the dipole interaction is mainly repulsive. 
However, when we decrease the value of cj^ the attrac- 
tive part of the dipolar interaction can give rise to pairing 
and superfluidity if temperature is sufficiently low 
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The radial quadrupole mode might then become a good 
test for the transition. For highly pancake traps the fre- 
quency should correspond to the coUisionless prediction 
and hence be given by the upper (red) line in Fig. [3] 
(bottom panel) . As the trap aspect ratio is decreased, at 
some point the attractive dipolar interactions will induce 
the gas to enter the superfluid phase, and the quadrupole 
frequency should drop to its hydrodynamic value, given 
by the lower line (black) in the figure. To give an order 
of magnitude, following Ref. [215 fo^' ^ Fermi gas contain- 
ing N — IQ^ atoms or molecules for which = 0.5, 
the transition to the superfluid phase should occur at 
Acr ^2: for A > 2 the gas would be in the normal phase, 
while for A < 4 it should be a BCS superfluid. 

Alternatively, in the normal phase at finite T, the 
study of the quadrupole oscillation and the transition 
from the coUisionless to the hydrodynamic regime could 
provide an important test of the role of the collisions 
caused by the dipolar interaction. Indeed, if r is the 
typical collision time due to the dipolar interaction and 
Lo the colective mode frequency one gets the coUisionless 
(hydrodynamics) result provided lot 3> 1 (wr ^1). If 
we estimate r at low temperature within the Bor n ap - 
proximation for the scattering amplitude (see e.g |23l). 
we get 



'dd' 



T 



-dd 



(6iVA)i/3 \huj 



(33) 



where, in the last equality we replace ep (and Tp) with 
its non-interacting value = hui±{&N\Y^^ . Since the 
frequency of the quadrupole oscillation is of order a;_L, 
the condition for being in the coUisionless regime is then 
equivalent to requiring T < Tp / {eddi'oN XY/*^). 



IV. VIRIAL THEOREM AND THE 
MEASUREMENT OF THE DIPOLAR 
INTERACTION 



In Sect. II the radial virial relationship Eq. ([29)) has 
been used to calculate the quadrupole mode frequency 
in the coUisionless regime. Such expression can be ob- 
tained by considering a radial scaling transformation of 
the form x — > e'^x, y — > e'^y applied to the ground state, 
with the proper normalization constraint (e.g., the den- 
sity changes as n — >■ e~^'^n(e'^a;, e'^y, z)) and by imposing 
that the total energy ([3]) vanishes at first order in c. A 
similar virial relation can be found by considering a scal- 
ing transformation along the z-th axis and one finds: 

2Ski„,. - 2^ho,. + ^dip + ^NepEddhin) = (34) 

By summing the radial and the longitudinal virial expres- 
sions Eqs. (P^j) and one recovers the more familiar 
virial expression 2i?kin — 2_Eho + 3£'dip = (see, e.g., 

ilH). 

The virial expressions play an important role also in 
the dynamics of the system following the sudden switch- 
ing off of the dipolar interaction. Let us suppose that 
the system, initially at equilibrium in the presence of the 
dipolar interaction, is perturbed at t = by suddenly 
setting Vdd = (in the case of heteronuclear molecules it 
is enough to set the external electric field equal to zero). 
For t > the dynamics of the system is governed by the 
single particle Hamiltonian Hq = /2m -f Vho{'^)- For 
example the average radial square radius will evolve in 
time according to the law 



m 
4 



Pi 
2m 



-{muj^rj_) 



(35) 



where, in the last identity, we have introduced the en- 
ergy E'l = {p]_/2m) + {mix>\r\l2) along the radial di- 
rection. This quantity is conserved in time (for t > 0), 
due to the absence of the dipolar interaction. After some 
straightforward algebra it is easy to find the following 
time dependence for the radial oscillator energy: 

-Bho,i(i) =i^ho,i(i = 0) 
+ [-Ekin,±(i = 0) - -Bho,i(t = 0)] sxx^iujxt) (36) 



Thanks to the virial theorem (|29l) the quantity in the 
square brackets in front of sin^(ajxt) can be directly ex- 
pressed in terms of the dipolar contribution calculated 
before switching off the dipolar force. A similar relation- 
ship holds also for the time evolution of the oscillator en- 
ergy along the z-th direction - same Eq. (j36]) with _L— > z 
- so that, combining the measurements of the time evolu- 
tion of the radii along the radial and axial directions (and 
hence of i?ho,±(i) and i?ho,z(i) respectively) one can even- 
tually determine experimentally the value of the dipole 
energy as well as the function /i(k) p3 |. 
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V. CONCLUSIONS 

In this paper we have studied the effect of the dipolar 
interaction on the surface quadrupole mode in a trapped 
Fermi gas in both the cohisionless and the (superfluid or 
coUisional) hydrodynamic regime. The predicted values 
depend on the strength of the interaction and the cloud 
aspect ratio. We have in particular found that: 

(i) in the hydrodynamic regime the frequency is af- 
fected weakly and is always smaller than the "standard" 
result \f2ijj}_ which is recovered for a strongly deformed 
cloud, i.e., when the anisotropic nature of the dipolc- 
dipole potential is irrelevant; 

(ii) in the cohisionless regime the frequency of the 
quadrupole mode can be strongly affected by the dipo- 
lar interaction and it is smaller (larger) than the free gas 
case 2wj^ for a pancake (cigar) shape trap. 

The difference between the frequencies in the two 
regimes is large enough to be experimentally resolved. 

It is important to notice that the gas can be in the 
hydrodynamic regime either because of superfluidity or 
because of collisions. We provided an estimate for the 
temperature in order the system be collisional for the 
qudrupole mode. 

The virial theorem, which has been employed to de- 
rive the quadrupole mode frequency in the cohisionless 
regime, has been also proved to be useful to measure the 
dipolar energy of the system via a sudden switch-off of 
the dipolar interaction itself. 
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where F{(p\a) and E{(p\a) are the incomplete elliptic 
integrals of first and second kind [2^ and their arguments 
are given by: 



sin If = ^1 - kI . (A4) 

Much information about this function can be found in 
Ref. [l^. It can be seen that when the two arguments 
are equal, = Ky, one recovers expression (|A1[) . Figure 
|4]shows the behavior of /(k) with respect to its argument. 
Note that it vanishes for k = 1, implying thus that the 
dipolar energy vanishes for a spherical dipolar gas. It is 
also interesting to notice that in the limit k — > oo the 
hmction tends to a constant value, /(k) — > —2, and the 
dipolar energy depends on the Thomas-Fermi radii in 
the same way as does a contact interaction. This is the 
underlying reason why the quadrupole frequency is not 
affected by dipolar interactions in the HD regime, when 
the gas approaches the 2D limit. 




FIG. 4: Anisotropy function /(k) and its derivatives, /i(k) 
and g{K), that appear in the expressions of the quadrupole 
frequency, as a function of k. 



Appendix A: Anisotropy functions /(k), g{ft) and /i(k) 

a. The function /{k) The function /(k) reflects the 
anisotropy of dipolar interactions both for bosonic and 
fcrmionic systems. For a cylindrically symmetric system 
it can be written as [201 



b. The function g{K) In the expression of the 
quadrupole frequency, in both regimes, there appears the 
function ^(k) defined by 



d^f{e^K,e-'>K) 



6=0 



1 + 2k^ 



I 



k2 (1 - 



tanh'^ - (Al) _ Sk^ ^1 - ^2(2 + ISk^) - k^{12 + 3k^) tanh"^ Vl - 



(1 - ^2)7/2 



If the system has no cylindrical symmetry, the anisotropy 
function depends on the aspect ratios = Rx/ Rz and 
Ky = Ry/Rz and has the form [l^ 

/(kx, i^y) = 1 + inxKy }" }- , ^ } , (A-2) 



(A5) 



It is shown in Fig. It vanishes for k — > and k — > oo 
and has a maximum at k ~ 2.37. 

c. The function h^n) In the radial and longitudinal 
virial expressions we have introduced the function /i(k) 
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defined by 



h{K) 



(1 - k2)5/2 



-3^1 - k2 + (2 + K^) tanh"^ Vl - 

(A6) 



The behavior of the function is very similar to the 
one of gin). It vanishes botli for k — > and k — >■ oo and 
has a maximum at k ~ 1.56 (see Fig. |4]). 
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